Maximally Tight Bell Inequalities Involving Lower Order Correlations 



Marcin Wiesniak,^ Mohamed Nawareg,^ and Marek Zukowski^ 

'institute of Theoretical Physics and Astrophysics, 
University of Gdansk, PL-80-952 Gdansk, Poland 

Most of known multipartite Bell inequalities are for correlation functions for all subsystems. We 
propose a method of derivation of tight Bell inequalities involving also correlations of lower order. An 
example of a A'' qubit state, which violates such inequalities despite vanishing A''-qubit correlations 
is given. We show how to derive a necessary condition for a state to violate such inequalities. 
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Entanglement might be seen as the most striking fea- 
ture of quantum mechanics With perfect correlations, 
a feature of maximally entangled states Einstein, Podol- 
sky, and Rosen 0| tried to argue that quantum theory 
could be supplemented with elements of reality respect- 
ing locality. This was disproved by Bell [3], who derived 
a theorem showing a discrepancy between all possible lo- 
cal realistic theories and entanglement. Since then Bell 
inequalities have been further generalized to more com- 
plicated experimental situations, and found applications 
in quantum cryptography, reduction of communication 
complexity, and other quantum information protocols. 
Series of new inequalities for an arbitrarily large collec- 
tion of qubits, with two alternative observables per side 
were derived in Refs. 4, 5, 7-9]. Basing on similar ideas 
inequalities with more settings per site were derived in 



' ^' Refs. [iSE 



Parallel to these developments, geometrical aspects of 



Bell inequalities were studied [1^. The set of all possi- 
ble probabilities of experimental results, for experiments 
involving a finite number of settings, describable by lo- 
cal realistic theories forms a convex "Pitowsky" polytope 
in a vector space. The dimension of the space is deter- 
mined by the number of particles, the dimensionality of 
the local observables, and by how many observables each 
observer is allowed to choose from. In case of iV-partite 
this will he D = 2^*=! "nn case of inequalities presented 
here, D = Y\^{ni 4- 1) — 1, with Ui being the number of 
observables at disposal of the ith observer. The extreme 
points, or vertices, of this polytope are given by "deter- 
ministic models" , that is apply to theoretical situations 
in which just one set of results is always obtained. For 
qubits, we have K — 2^^=!"'. A tight Bell inequality 
is a one linked with a hyperplane spanned by at least D 
vertices of the polytope, all lying on the same facet, and 
thus defining it. All such vertices saturate the inequality. 
Here we would like to introduce a notion of maximally 
tight inequalities. We shall use this term whenever the 
associated hyperplane contains half of the vertices of the 
polytope. An example of such an inequality is the cel- 
ebrated CHSH inequality [13[, which takes value +2 for 
eight deterministic local realistic models, and —2 for the 
other eight. Basing on such concepts one can develop 
methods to derive maximally tight Bell inequalities for n 
qubits with three observables per site 1^ . Some explicit 



forms of such inequalities were obtained in Ref. |15| . 

One can find multi-qubit states, which, despite be- 
ing entangled, do not posses any correlations between all 
qubits [16]. In fact, they have non-zero correlations only 
between an even number of particles. Such states are 
now within the experimental reach, see for Dicke states 
refs. [l3j[3| or for rotationaUy invariant states refs. [I^. 
Naturally, as there are no correlations between all qubits, 
such states cannot violate any of the aforementioned in- 
equalities, as they involve correlations of the highest or- 
der. 

Here we present an algorithm to derive tight Bell in- 
equalities. These inequalities involve only some correla- 
tion functions for all N qubits, and mostly rely on lower 
order (in particular, [N — 1) qubit) correlations. We will 
demonstrate that they can be violated by states without 
any A^-fold correlations, for an odd N . 

Let us start with the simplest case, for which one can 
derive such inequalities, that is for three qubits. We as- 
sume that each observer can make one of two dichotomic 
measurements. The local realistic values, for the deter- 
ministic case, will be denoted as A\, where i = 1,2,3 
denotes the observer and fc = 1, 2 is an index of the local 
observable. We put, as usual, the possible actual hidden 
values as A\. = ±1. The set of inequalities that we shall 
derive is based on the following algebraic identity: 

^iAl + Al){Al + Al)iAl + Al) 
±1{AI + Al){Al - Al) ±\{Al + Al){Al - A^) 



±\{A\-Al){Al 



±UAl-A^,){Al-Al){Al-Al)^±l, 



(1) 



which we shall prove below, not by inspection, but using 
the general method that we want to present here. 

The left hand side of (H]) can be expressed as a scalar 
product of two tensors. The first one is built out of (de- 
terministic) local realistic predictions, and reads 

HA) = {AlAi 1) ® (Al All)® (AlAl 1), (2) 

where A^ {Al, A^, Af, A^, Af, A^). Note that there are 
2^ = 64 tensors of such a kind. They are in a one-to-one 
way relation with the vertices of the Pitowsky convex 
polytope representing all possible local hidden variable 
representations of the problem, and thus we shall call 
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them simply vertices. The tensor component which is 
always equal to 1 plays here only a technical role. 

Consider now another tensor built out of linear com- 
binations of the following product tensors: Vbib2b3 = 
Vbi ®Vb2 ®Vb3, where bi = +1, —1, 0, and v±i = (1, ±1, 0) 
and fo = (0,0, 1). To simplify notation we shall denote 
the first two vectors as v±, an use a similar convention 
for the tensors Vbib^bs- Whenever we restrict the values 
of 6j;'s to ±1, or ±, we shall denote them as a^. The 
vectors form a unnormalized basis in R^. Next we con- 
struct a specific linear combination of Vbib2b3 which has 
the property that its scalar product with any vertex has 
the same modulus, in our case 1, and additionally only 
one term of the linear combination has a non-zero scalar 
product with any of vertices. To this end, one must mul- 
tiply any Wai02 03 by |. Next we observe that such tensors 
have 8 = 2^ vertices with which they are non-orthogonal. 
Those w's with just one bi — we multiply by j. Each 
of them has two times more vertices, than the tensors of 
the earlier class, which are not orthogonal. And so on. 
Of course, vqqq cannot enter a non-trivial construction of 
such a kind, as they are non-orthogonal with respect to 
all vertices. 

Our final task is to find a set of w's having the required 
property. We can start with and arbitrary one of the kind 
^'±±±1 a-nd if we continue only with such ones, we shall 
end up with a WWWZB inequality. To get anything new 
we have to use at least one v with at least one index. 
So let us start with "0+++, and in the second step add 

Q, as they certainly do not share any non-orthogonal 

vertex (generally Vbib2ba f^nd "S^^b^ do not share a non- 
orthogonal vertex, if they have a pair of indices satisfying 
^ ¥^ bk b'f. ^ 0) . A similar situation happens if we add 

in the second step woh or V-o+. Next notice that all of 

them do not share non-orthogonal vertices, so the whole 
set is in. Finally we see that we could have started with 

V to get the same triad with zeros, and that this 

vector does not share any non-orthogonal h with "0+++. 
As there are only 64 vertices and 2x8 + 3x16 = 64, we 
have built a linear combination of the required properties, 
which reads 

± \(^+++ =^ ^ ) =^ ^("+-0 ± ^-0+ ± '^0+-)- (3) 

Obviously a "sum of products of components" scalar 
product of © and ^ gives (P). 

Note important facts. The coefficients give us fraction 
of vertices of the Pitowsky polytope with which the given 
member of the combination is non-orthogonal (e.g., t)+++ 
is non-orthogonal to i of the vertices). Note further, 
that as the moduli of the coefficients add up to 1 by 
putting 'C'(&i62&3) = ±i'(&i62&3), where ± now stands 
for the effective sign in the combination before the given 
tensor, we have a convex combination of the form 

^(w+++ + v'—) + |(w+-o + ^^-0+ + ^^0+-)- 



After the usual averaging we get a set of 2^ Bell in- 
equalities (inequalities of this kind were found in ref. (20j 
using a completely different method, not allowing for easy 
generalizations) 

±{^{A\ + Al)iAl + Al){Al + Al) 
±1{A\ + ADiAi - Al) ± liAi + ADiAl - A^) 
±l{Al-Al){Al + Al) 
±l{Al-Al)iAl-Al){Al-Al))<l, (4) 

which just as in the case of WWWZB ones can be 
wrapped up to just one (this feature was not noticed in 

\{l{A\+Al){Al + Al){Al + Al))\ 

+|(i(A} + Ai)(Af - Al))\ + \{\{Al + A^,)iAl - Al))\ 

+ mAl-Al){Al + Al))\ 

+ \{^{Al-Al){Al^Al){Al-Am<h (5) 

Unfortunately ineq. ([5]) is not violated by any state. Still 
by extending the above method to more qubits we get 
"relevant" tight inequalities. Note further that by a con- 
sistent sign swap ± to =F for one particle, or more of them, 
in all terms of ([5]), we get other inequalities. 

Let us now describe a general method for the N qubit 
case. We define Z(N, k) as a certain set of A^-fold ten- 
sor products of u_ and wq involving at most k wq's, 
We shall again construct an algebraic identity, which can 
be interpreted as scalar product of a tensor which is a 
linear combination (interpretable as a convex one) of the 
elements of Z{N, k), with another tensor, built out of lo- 
cal realistic values, of the form h{A) — (E)fLi{A\, A2, 1), 
where now A = {A\, A^, A^ , A2 )■ Obviously, the lat- 
ter tensor is in a one-to-one relation with a vertex of the 
Pitowsky polytope of the problem. We shall call it vertex. 

Although this is not necessary, we focus on Z{N,kys 
which are invariant with respect of a cyclic permutation 
of the indices of the tensors. We shall call such set as 
having a cyclic permutation symmetry (CPS). 

The set Z{N, k) is constructed in an iterative way. We 
can start to build Z{N, k) with an initial set of tensors of 
the form Vaia2... — Vai ®Wa2---j hereafter denoted as w±±..., 
each of which has a CPS. We extend the definitions given 
in the example to N particles, thus as before — ±1, 
etc. For prime values of iV, there are only two such vec- 
tors and V Alternatively we could choose 

one of the proper subsets with CPS. In next steps, we 
supplement Z{N, k) with specially chosen tensors Vbi...bN 
with at least one and at most k of biS equal 0. Before we 
proceed further we notice the following: vertices can be 
grouped into classes of equivalence. The following prop- 
erty the scalar products of a tensor Vaia2,...aN with the 
vertices define such classes: there are groups of 2^ ver- 
tices which have a non zero scalar products with it, equal 
to ±2''^. All other vertices are orthogonal to it. Such a 
class will be denoted by ft(a), where a = (ai, ..ajv). 
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We repeat the following steps: 

1. We make a table of scalar products between the 
vertices and vectors belonging to the current set 
Z{N, k). 

2. Out of all groups of vertices h{a) orthogonal to all 
V in the current Z{N, k) we choose 2'^ such, which 
are different from each other at not more than k 
fixed positions (in the earlier example this would 

have been e.g. h{-\ h) and h{-\ )). We take 

associated with them 'Caia2,...ajv's and replace u+'s 
and V-S with vqS at these chosen positions (in the 
example, we would obtain v+_o). 

3. We add to Z{N, k) all cyclic permutations of the 
tensor obtained in the last step (involving Wq's). 

This procedure is repeated until Z{N, k) is complete, i. 

e., there are no vertices orthogonal to all elements of 
Z{N,k). The structure of Z{N,k) must be such that 
each vertex has a non-zero scalar product with exactly 
one tensor v G Z{N, k). Only then it generates a maxi- 
mally tight inequality. It may also happen that there are 
not enough vectors to complete step 2. This can hap- 
pen, for example, while constructing Z(7, 3). One could 
then temporarily decrease k and continue the protocol, 
or simply leave Z{N, k) incomplete, at the price of not 
getting a maximally tight inequality . 
Finally we build a tensor: 

Z{Z{N,k))= S{b)2-^^'>v{b), (6) 

Z{N,k) 

where B{b) = J2i=i N and S{b) is any sign function 
of the indices, that is of two values ±1. Note that the 
factor 2~^^''^ guarantees that the term containing v{b) 



For five qubits we have found 



contributes to the scalar product with any vertex ±1 or 0. 
The properties of the set Z{N, k), provided it is complete, 
in turn guarantee that for any vertex h{A) 

Z{Z{N,k))-h{A)=±l, (7) 

where the scalar product, denoted a •, is a "sum of prod- 
ucts of components" one. This after averaging over any 
possible probabilistic distributions of vertices, gives a set 
of 2^ inequalities, where z is the number of elements 
in Z(N, k) (this reflects the number of possible relative 
signs between the elements of the tensor Z). One can 
wrap them into a single inequality using exactly the same 
method as in the case of the example for A'' = 3 given 
earlier. 

Some remarks on the construction. The first observa- 
tion is that, since the protocol does not determine which 
vectors will be chosen in each step, but rather leaves it 
up to us, there might exist different sets Z{N, k) for the 
same values of N and k leaving to different inequalities. 

The second point is about values of k. Certainly, we 
are interested in odd fc's and odd A'''s. There are states 
which have correlations of even orders, but with no of 
odd ones. States with the opposite property are unknown 
(probably impossible). One must have A: < ^. This is 
because otherwise, there could be some vertices, which 
have a non-zero scalar product with more than one vector 
belonging to Z{N, k). Too see this, e. g., take a pair of 
t;+oo and vo()+ of the N = i example. In such cases the 
construction cannot work. 

Let us now present some specific sets Z{N, k) and dis- 
cuss their violation (CP denotes the set of all cyclic per- 
mutations, e.g., CP{v^ o) = {^'^ Oi^oh ,w_o+})- As it 

has been said above the interesting inequalities are those 
for odd N. 



Z(5, 1) = {v+++++,v } U CP{v+++.o) U CP(z;_++_o) U CP(i;_+__o), (8) 



for seven 



Z{7, 1) = V } U CP(i;+++++_o) U CP(«+++___o) U CP(.;_+__++o) U CP{v +o) 

UCP{v+-+.++o) U CP(i;++_++_o) U CP{v.-+-++o) U CP{v+.+ o) U CP{v+.++.+o), (9) 
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and for nine, which is a non-prime number 

Z(9,l) = {«+++++++++, w 

■i^ + + +,V++- + + - + + -,V+-++-++-+,V-++-++-++} 

U CPiv+++++++-o) U CP(«_++++++_o) U CP(«_+++++_+o) U CP(«_++++_++o) 

U CP(i;_+++++__o) U CP(«_++++_+_o) U CP(v_++++__+o) U CP(«_++++___o) 

-+++-H a) 

U CP(w_+++__++o) U CP(w_+++ o) U CP(i;_+++___+o) U CP(i;_++_++__o) 

U CP(«_++_+_+_o) U CP(i;_++_+___o) U CP(«_++__++_o) U CP(i;_++ o) 

U CP(«_++___+_o) U CP(«_++_+__+o) U CP(«__++__+_o) U CP(«_+_+__+_o) 

U CP(«_+_+ o)UCP(i;__+__+_+o)UCP(z; U CP(«___+ „), (10) 

Z(9,3) = {«+++++++++, u ,u+ — + — + — ,v-+ — + — 

V + + + , + + - + + + + + + 

U CP(v+++++_ooo) U CP(w_++++_ooo) U C7P(v_+++__ooo) U CP(w_++___ooo) 
U CP(?;_++o-+-oo)UCP(w_++_+_ooo)UCP(i;__+_o-o-o)- (H) 

I 



We have studied violations of these inequalities, as- 
suming a symmetry between observers. We took as lo- 
cal observables A^^ = cos (piUz + sin ^iCTj., and A2 = 
cos (piaz —sin 4>iax (for an explanation of such a choice, see 
further), with all 0i = 0. The maximal eigenvalue of for 
five-qubit Bell operators related to Z(5, 1) was found to 
be 1.97435 with 0i = f . For Z(7, 1), Z(9, 1) and Z(9, 3) 



we assumed (pi = (j) and the maximal observed violation 
ratios were 1.84331, 2.18414 and 1.79497, respectively. 
One of the five-qubit states providing the maximal vi- 
olation of the wrapping inequality with moduli defined 
by Z(5, 1), for Pauli observables constrained to the XZ 
plane, however without the (f>i = (j) constraint, reads 



lip) « a|00000) + 6(-|00011) + |00110) + |01100) - |11000) - |10001)) 
-l-c(lOOlOl) - 101010) + llOlOO) - 101001) - 110010)) + d(|00111) -t- |01110) + 111100) - 111001) - 110011)) 
-t-e(lOlOll) - 110110) - 101101) + IllOlO) - 110101)) + /(|01111) + |10111) - |11011) + |11101) + |11110)) 

+5IIIIII), (12) 



where a = 0.462854,6 = 0.161096, c = 0.19409, d = 
0.181191, e = 0.220891, / = 0.107669, and g = 0.039699. 
The state is given in the simplest form, a generalized 
Schmidt decomposition, defined in 22 1. 

We can apply to this state a NOT map, which flips the 

NOT 

Bloch vector of a qubit state (CTa — > — (Ta, a = x,y,z), to 
all qubits. Such a map is composed of unitary rotations 
and a full transposition, and thus it is positive one. If we 
take an equal mixture p = i(|V')(V'| -t- NOT'»^{\ip){'ip\)), 
the corresponding mean value reads 1.798, still much 
above 1. Notice that p possesses the same 4-party corre- 
lations as |'0)(V'|i but it is globally uncorrelated. Hence 
it could not be used to violate any non-trivial WWWZB 
inequality for five qubits . 

Let us now briefly present a method to establish a nec- 
essary condition to violate the inequalities. We shall use 
the approach of Ref. In the course of constructing 



the inequalities the pairs of local observables occur ei- 
ther as ^{A\ -I- A^) or ^{A\ — A2). Take any possible 
pair of observables for the party i, that is A'l = n\ ■ a 
and A2 = n2 ■ (T, where n\ are unit vectors, and (? is a 
quasi- vector formed of Pauli operators for i. As n\ + 
and n\ — n\ form an orthogonal pair, after some alge- 
bra we see that we can put ^(A^^ -I- A2) — cos piw"^ ■ a 
and ^(A^^ — Aj) = sin^iW^' ■ ct, where pi is certain an- 
gle, dependent on n\ and and and w^"^ form an 
orthonormal pair. For simplicity we put = and 
_ TiO^ can rewrite our Bell expression (for 

technical details of the approach consult 9]) for the wrap- 
ping inequality with moduli as a scalar product between 
two vectors, one containing trigonometric functions, and 
the other moduli of elements of the correlation tensor, 
|Py...| = \Trp{(Ji ® ...)|, with the indices extending 
from to 3, and ctq = I. The inequality [5] for three qubits 
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(which as aforementioned cannot be violated) gives 

( I cos 6*1 cos 62 cos 6*3 1 , 1 sin 9i sin 62 sin 6*3 1 , 
I cos6'i sin02|, I COS02 sin^sj, | sin 6*1 costal) 



cus ui sill 1/2 I J I i^'-'s 'J2 f "3 1 1 I '-'1 '^3 
■(|r333UTin|,|r310|,|T03l|,|ri03|) < 1. 



(13) 



It is now easy to see that the norm of the first vector is 
at most 1 (the equahty holds for the complete Z{N, k)). 
Thus, Cauchy-Schwartz inequality gives 



/t-t2 

^333 



^310 



^031 



7^103 <1: 



(14) 



as a sufRcient condition to satisfy the inequality. If it 
holds for a given three-qubit state for all local bases, no 
inequality generated by this Z can be violated by this 
state. 

Similarly, the sufficient condition for a state to satisfy 
all inequalities generated by ,e.g., Z(5, 1) is given by 



+ 



^33333 
q-,2 
33103 



^niii 

^W310 

rp2 I rp2 

^13110 + -'01311 



^33310 
^01331 



^03331 
^10133 



^n)333 + ^31033 



^31013 



33101 



rp2 

10131 



^noi3 + ^3^101 ^ 1- (15) 



The presented method leads to tight mutlipartite mul- 
tisetting Bell inequalities which are an extension of the 
WWWZB ones to multi order correlations. The actual 



construction was inspired by Ref. 15|. Our results cer- 
tainly opens many interesting questions requiring further 
investigation. For example, how many different Z's can 
be constructed for given iV, or what other symmetries 
can be imposed? A different task is construction with 
this method of inequalities violated by measurement data 
for experimentally accessible states. The degree of rele- 
vance of these inequalities in comparison to the known 
ones is another class of problems worth investigating. 

While preparing this work we have learned about a new 
class of multisetting Bell inequalities for three qubits '21'] . 
It is certainly interesting how these inequalities generalize 
to more qubits and if the two families, the one emerging 
from Vertesi-Pal inequalities and the one derived from 
here, converge. 

This work is part of EU program QESSENCE (Con- 
tract No.248095) and MNiSW (NCN) Grant no. N202 
208538. MN is supported by the International PhD 
Project "Physics of future quantum-based information 



technologies" : grant MPD/2009-3/4 from Foundation for 
Polish Science. 
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